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ABSTRACT 


Like  the  weather,  upper  atmospheric  density  often  behaves 
unpredlctably.  Even  the  best  models  of  the  upper  atmosphere  cannot 
consistently  predict  atmospheric  density  for  Transit  System  satel¬ 
lites  (1000  km)  to  the  desired  high  precision. 

A  comprehensive  evaluation  Is  presented  of  the  drag  estima¬ 
tion  concept  as  a  tool  In  routinely  determining  the  orbit  of  the 
Transit  System;  the  benefits  to  the  Transit  System  are  also  exam¬ 
ined.  The  resulting  data  (fitted  along-track  bias  force)  can  be 
used  subsequently  to  study  upper  atmospheric  density  variations  at 
Transit  satellite  altitudes  that  will  lead  ultimately  to  Improved 
upper  atmospheric  density  models. 

Drag  estimation  In  satellite  orbit  determination  yields 
higher  precision-tracked  (determined)  ephemerides  and  significantly 
Improves  the  prediction  (extrapolated  satellite  position)  accuracy 
of  the  Transit  System. 
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1.  INTRODUCTION 


It  has  become  apparent  that  a  complete  and  accurate  repre¬ 
sentation  of  the  effects  of  drag  on  satellites  is  lacking.  Even 
the  best  models  of  the  atmosphere  cannot  duplicate  the  actual  pre¬ 
vailing  conditions  at  satellite  altitudes.  An  alternative  (or 
addition)  to  modeling  the  atmosphere  Is  to  eliminate  Its  effects  on 
satellites  by  using  "drag  estimation,"  which  Is  the  inclusion  of  a 
drag  parameter  in  the  satellite  orbit  determination  process.  The 
resulting  fitted  parameter  Is  used  to  generate  an  along-track  force 
correction,  thereby  compensating  for  the  shortcomings  of  the  drag 
model. 

Errors  In  the  modeled  drag  force  result  In  along-track  satel¬ 
lite  position  errors  that  grow  quadratlcally  with  time  (see  Appen¬ 
dix  A).  Figure  1  Illustrates  some  typical  fixed-site,  along-track 
navigation  results  that  are  proportional  (and  nearly  Identical)  to 
the  satellite  along-track  position  errors.  The  orbit  determination 
Interval  gives  a  measure  of  the  satellite  position  errors  after  an 
orbit  adjustment.  The  ephemeris  extrapolation  interval  shows  a 
quadratic  growth  In  time  (centered  at  the  middle  of  the  track)  and 
gives  a  measure  of  satellite  position  errors  before  an  orbit 
adjustment. 

Without  drag  estimation,  the  theoretical  function  used  in 
the  orbit  determination  process  includes  constant  and  linear  (In 
time)  terms.  The  resulting  fitted  parameters  are  converted  to  cor¬ 
rections  to  the  satellite  semimajor  axis  (linear  term)  and  mean 
anomaly  (constant  term).  Nothing  has  been  done  about  the  quadratic 
growth  that  is  evident  In  the  track  and  that  results  in  serious 
errors  at  the  end  of  the  prediction  span.  Drag  estimation  includes 
In  the  fitting  process  an  additional  term  that  grows  quadratlcally 
In  time.  The  resulting  fitted  parameter  can  then  be  converted  Into 
an  along-track  force  correction  (correction  to  the  modeled  drag 
force),  which  is  handled  like  any  other  force  acting  on  the  satellite. 


residuals  (m) 
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Fig.  1  Along-track  navigation  residuals  for  satellite  30120  in  1978. 
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2.  SUMMARY  AND  CONCLUSION 


Drag  estimation  is  a  viable  concept  that  can  be  implemented 
and  used  routinely  in  the  orbit  determination  process  for  the  Tran¬ 
sit  System.  The  primary  assumption  that  forms  the  basis  for  the 
drag  estimation  technique  is  that  drag  model  errors  vary  slowly  with 
time;  that  is,  correlation  time  is  long  compared  with  the  orbit  de¬ 
termination  Interval.  Results  presented  here  confirm  that,  on  the 
average,  the  assumption  is  valid.  Occasionally,  however,  rapid 
(compared  with  the  orbit  determination  interval)  changes  in  the  at¬ 
mosphere  result  in  poor  ephemeris  extrapolation  unless  the  result¬ 
ing  estimated  along-track  correction  force  is  intentionally  damped 
to  suppress  the  short-term  changes  In  the  atmosphere. 

A  manual  drag  estimation  system  has  been  used  successfully 
at  the  Navy  Astronautics  Group  (NAVASTROGRU)  since  April  1978.  The 
manual  procedure  involves  the  routine  monitoring  of  satellite  along- 
track  errors  (examination  of  predicted  fixed-site  along-track  navi¬ 
gation  errors)  and,  when  necessary,  the  correction  of  drag  force  by 
changing  the  coefficient  of  drag  (C^)  from  its  "true"  value  of  3.0. 

Automated  (computer-generated)  drag  estimation  removes  the  guesswork 
(magnitude  of  needed  changes  to  C^) ,  follows  short-term  (dally)  as 

well  as  long-term  drag  model  errors,  and  requires  substantially  less 
manual  interference.  Benefits  include  higher  precis ion- tracked  eph- 
emerides  (sometimes  referred  to  as  precise  ephemerides) ,  Improved 
prediction  accuracy  (when  using  the  current  prediction  span) ,  and 
lengthening  the  prediction  span  without  loss  in  (current)  system 
performance.  Routine  drag  tracking  produces  valuable  data  (fitted 
bias  forces  and  mean  modeled  densities)  that  are  essential  for  the 
study  of  upper  atmospheric  density  variations. 
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3.  BACKGROUND 


Drag  estimation  was  first  attempted  In  1968  In  a  special 
version  of  the  Transit  System's  Orbit  Improvement  Program  (OIP) . 

Its  purpose  was  to  supplement  the  existing  drag  model  and  thereby 
improve  orbit  prediction  accuracy.  The  drag  parameter  chosen  was 
C^,  which  for  the  Transit  satellites  had  a  nominal  value  of  3.0. 

Modeled  density  errors  would  be  countered  by  adjusting  (fitting) 
either  upward  (model  underestimating  density)  or  downward  (model 

overestimating  density)  from  its  nominal  value.  The  necessary 
theory  (see  Appendix  A)  was  developed ,  and  a  special  drag  tracking 
version  of  OIP  was  generated  and  tested  shortly  thereafter  on 
NAVASTROGRU's  IBM  7094  computer. 

The  conclusions  drawn  from  that  first  study  are  valid  today: 

1.  Drag  estimation  requires  data  intervals  longer  than  the 
currently  used  36  h  spans. 

2.  Even  with  longer  spans  (44  h),  the  variation  from  one 
span  to  the  next  of  the  fitted  values  of  was  too  sporadic  to  be 

attributed  to  the  uncertainty  in  the  drag  model. 

3.  The  linear  and  quadratic  coefficients  in  the  fitted  func¬ 
tion  (a^  and  a^)  cannot  be  treated  Independently  as  was  believed 

previously. 

Dramatic  improvements  in  the  tracking  accuracy  of  the  Transit 
System  were  achieved  during  the  intervening  years,  but  repeated  at¬ 
tempts  to  fit  a  drag  parameter  have  failed  to  alter  the  1968  con¬ 
clusions.  The  crux  of  the  problem  was  the  inseparability  of  the 
linear  and  quadratic  terms  (in  time)  in  the  orbit-fitting  equation. 
The  concept  of  fitting  out  a  drag  parameter  was  later  replaced  with 
a  bias  force  in  HLC  coordinates  (Appendix  A)  whose  along-track  com¬ 
ponent  (L)  can  be  adjusted  (on  option)  along  with  the  other  fitting 
parameters.  The  same  approach  was  also  carried  over  into  the  new 
Orbit  Determination  Program  (ODP),  the  IBM  360  replacement  for  OIP, 
where  the  ability  to  fit  out  an  along-track  bias  force  was  retained. 
Recent  increases  in  the  sun's  activity  and  the  associated  variations 
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In  upper  atmospheric  density  have  again  surfaced  to  remind  us  that 
the  upper  atmosphere  Is  unpredictable  (Ref.  1). 

This  report  reexamines  the  old  approach,  Introduces  a  new 
one  and  summarizes  the  results  of  extensive  testing  of  drag 
estimation  at  APL  and  NAVASTROGRU  where  the  concept  Is  being  evalu¬ 
ated  for  Inclusion  in  the  routine  operations  of  the  Transit  System. 


Ref.  1.  H.  E.  Hlnteregger,  "Development  of  Solar  Cycle  21 
Observed  In  EUV  Spectrum  and  Atmospheric  Absorption,"  J.  Geophys. 
Res. ,  Vol.  84,  No.  A5,  1  May  1979. 
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4.  THE  LATEST  ATTEMPTS  USING  THE  OLD  APPROACH 


In  theory,  there  is  no  problem  in  separating  the  linear  from 
the  quadratic  growth  in  the  satellite  along-track  errors.  In  prac¬ 
tice,  we  have  to  contend  with  a  number  of  problems: 

1 .  The  data  are  noisy, 

2.  The  data  are  sparse  (small  number  of  passes)  and  poorly 
distributed  geographically,  and 

3.  Residual  pole  position  errors  introduce  a  24  h  oscilla¬ 
tion  into  the  along-track  navigation  residuals. 

The  approach  that  had  been  taken  previously  was  to  assign 
quadratlcally  growing,  along-track  position  errors  to  drag  model 
errors.  A  bias  force  in  the  along-track  direction  was  introduced 
to  compensate  for  the  modeled  errors.  One  difficulty  is  that  the 
quadratic  and  linear  terms  "fight"  each  other;  that  is,  it  is  quite 
common  to  fit  out  large  compensating  quadratic  and  linear  terms. 

The  linear  term  is  converted  to  an  adjustment  to  the  initial  semi- 
major  axis  (satellite  period),  and  the  quadratic  term  is  trans¬ 
formed  into  an  along-track  bias  force.  If  we  fit  out  an  incorrect 
(unusually  large)  quadratic  term  that  is  compensated  for  by  an 
equally  large  (and  opposite  in  sign)  linear  term,  we  will  pay  the 
price  for  these  errors  in  the  update  (ephemeris  extrapolation)  in¬ 
terval.  The  "incorrect”  bias  force  will  cause  quadratic  (in  time) 
growth  in  the  along-track  errors  that  will  only  be  partially  com¬ 
pensated  for  by  the  equally  incorrect  though  only  linearly  growing 
(in  time)  period  error.  The  result  could  be  disastrous.  In  a  re¬ 
cent  test,  two  out  of  three  consecutive  orbit  determination  runs 
resulted  in  poorer  predictions  when  compared  to  their  standard  op¬ 
erational  counterparts,  upf.ng  four  clusters  (48  h  of  data)  for  orbit 
determination  and  ephemeris  extrapolation  intervals  (see  Table  1). 

In  a  second  test,  we  added  two  more  clusters  (24  h  of  data) 
to  the  orbit  determination  intervals  but  kept  the  extrapolation 
spans  to  four  clusters  (48  h).  The  results  were  better;  two  out  of 
three  improved,  with  the  third  showing  only  a  slight  degradation 
when  compared  to  the  no-drag  fitting  case  (Table  2). 

In  the  runs  shown  in  Tables  1  and  2,  we  intentionally  intro¬ 
duced  an  error  in  the  drag  force  by  setting  the  mean  solar  index  to 
125,  which  Is  30  units  too  high.  Therefore,  we  would  expect  to  fit 
out  a  bias  force  in  the  positive  along-track  direction  (opposite  to 
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Table  1 

Drag  estimation  (includes  fit  to  linear  terms)  results  using  2  days  of  data. 


Satellite 


30140 


Span 

(1977) 

From 

To 

(day)  (h) 

(day) 

(h) 

339  18 

341 

12 

341  18 

343 

12 

343  16 

345 

12 

345  16 

347 

12 

Along-track 
bias  force* 

(R./s2  x  1015) 


rms  along-track  errors 
(m),  48  h  extrapolation 

Drag  fit 

No-drag  fit 

9.7 

29.7 

26.2 

15.2 

56.0 

28.7 

*RQ  scaling  factor  -  6378.166  km 


Table  2 

Drag  estimation  (includes  fit  to  linear  terms)  results  using  3  days  of  data. 


Span*  (1977) 
From  To 

Satellite  (day)  (h)  (day)  (h) 


30140  338  18  341  12 

340  18  343  12 

342  18  345  12 

344  16  347  12 


Along-track 
bias  force 

(R„/s2  x  io15) 


rms  along-track  errors 
(m),  48  h  extrapolation 

Drag  fit 

No-drag  fit 

12.8 

29.7 

13.6 

15.2 

35.1 

28.7 

*Note  that  the  extrapolation  span  covers  only  the  last  48  h  of  each  orbit 
determination. 
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th«  drag  force)  to  compensate  for  the  error  In  the  drag  force  (nega¬ 
tive  along-track  direction) .  It  was  surprising  to  find  that  In  the 
fourth  track  there  was  a  reversal  In  the  sign  of  the  fitted  bias 
force.  A  possible  explanation  will  be  discussed  In  the  last  section 
of  this  report. 


i 
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5.  A  NEW  APPROACH  TO  THE  OLD  PROBLEM 


The  difficulty  with  getting  drag  estimation  to  work  Is  the 
poor  (or  Incomplete)  separation  between  the  linear  (period)  and 
quadratic  (bias  force)  parameters  in  the  orbit  adjustment  process. 
The  linear  term  (period  error)  has  alway  been  one  of  the  standard 
six  orbit  parameters.  The  error  in  the  period  is  corrected  by  ad¬ 
justing  the  semimajor  axis.  Why  is  there  a  period  error?  Primar¬ 
ily  it  is  because  the  modeled  drag  force  is  less  than  perfect,  and 
the  period  (semimajor  axis)  needs  to  be  corrected  along  with  the 
other  parameters.  If  the  main  reasons  for  adjusting  the  semimajor 
axis  are  errors  in  the  modeled  drag  force  and  if  errors  in  the 
modeled  drag  force  result  in  quadratic  (in  time)  growth  of  the 
along-track  position  of  the  satellite,  would  it  not  make  more  sense 
to  fit  out  a  quadratic  term  and  forget  the  linear  term? 

To  test  the  validity  of  the  new  approach,  we  repeated  the 
four  orbit  determination  runs  and  compared  the  results  to  the  pre¬ 
vious  set  of  runs  (Tables  1  and  2) .  The  ODP  program  was  designed 
with  the  capability  to  determine  all  or  any  subset  of  11  parameters. 
Normally  (operationally)  we  adjust  eight  parameters  (six  orbit  and 
two  frequency  parameters) .  One  of  the  eight  (a^)  is  a  correction 

to  the  semimajor  axis  (or  the  linear  term).  With  our  new  approach, 
using  ODP,  we  remove  a^  from  the  fitting  space  and  include  instead 

the  along-track  bias  force  parameter  (or  the  quadratic  term,  a^). 

The  resulting  bias  force  is  then  routinely  integrated  along  with 
the  other  forces  acting  on  the  satellite.  Explicit  corrections  to 
the  semimajor  axis  (a^)  are  neither  computed  nor  applied.  The  re¬ 
sults  of  the  first  four  orbit  determination  runs  were  encouraging; 
two  of  the  three  extrapolated  ephemerldes  showed  significant  im¬ 
provement  over  both  of  the  previous  drag  estimation  attempts  (Tables 
1  and  2)  and  the  parallel  (OIP)  no-drag-fit  runs.  The  third  was 
only  slightly  worse  (rms  *  18.3  m  compared  to  15.2  m  (see  Table  3)). 

The  results,  which  warranted  the  continuation  of  the  experi¬ 
ment  for  an  additional  nine  intervals  (48  h  each)  or  a  total  of  26 
days,  were  better  than  we  had  anticipated.  Only  one  out  of  the  12 
extrapolations  was  significantly  worse  (30  m  rms  compared  to  15  m 
for  the  OIP  run) .  Eight  of  the  12  showed  significant  improvement 
(5  to  40  m  reduction  in  the  predicted  rms).  The  remaining  three 
did  about  the  same  (two  doing  slightly  better  and  one  slightly 
worse).  Table  4  summarizes  the  results.  A  few  observations 
follow: 


THE  JOHNS  HOSKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

LMJMCL  MARVIANO 


Table  3 

Drag  estimation  (excludes  fit  to  linear  terms)  results  using  2  days  of  data. 


Satellite 

Span 

(1977) 

Along- track 
bias  force 

(RQ/a2  x  1015) 

rms  along-track  errors 
(m) ,  48  h  extrapolation 

From 

(day)  (h) 

To 

(day) 

(h) 

Drag  fit 

No-drag  fit 

30140 

D 

B 

0.11 

— 

— 

H 

1 

0.15 

7.5 

29.7 

343 

16 

B 

0.03 

18.3 

15.2 

345 

16 

B 

-0.11 

16.2 

28.7 

Table  4 

Results  of  drag  estimation  (excludes  fit  to  linear  terms)  (ODP)  compared 
with  no-drag  estimation  (OIP). 


Satellite 

Span  (1977) 

Number  of  paaaea 

Fitted 

along-track 

force 

(R0/«2  «  1015) 
ODP 

48  h  along-track  residuals  (m) 

Orbit  determination  ephemerls 
extrapolation 

rms 

Last 

pass 

rms 

From 

(day)  (h) 

To 

(day) 

(h) 

Input 

Used 

(ODP/OIP) 

ODP 

OIP 

ODP 

OIP 

ODP 

OIP 

30140 

339 

18 

341 

12 

20 

20/20 

♦0.11 

3.9 

5.3 

— 

— 

— 

— 

30140 

341 

18 

343 

12 

25 

25/25 

+0.16 

4.2 

3.9 

+12.3 

+46.1 

7.5 

29.7 

30140 

343 

16 

345 

12 

27 

23/27 

+0.03 

3.0 

-31.2 

+23.1 

18.3 

B 

30140 

345 

16 

347 

12 

26 

26/26 

+0.11 

3.1 

3.2 

-30.7 

-54.6 

16.2 

E9 

30140 

347 

16 

349 

12 

25 

21/25 

+0.12 

3.4 

5.1 

+53.8 

-20.5 

30.6 

15.1 

30140 

349 

16 

351 

12 

24 

23/24 

+0.09 

3.2 

4.2 

+  8.0 

+26.1 

12.2 

16.1 

30140 

351 

16 

353 

12 

27 

25/24 

+0.08 

3.0 

3.9 

+16.8 

+28.2 

14.4 

19.4 

30140 

353 

15 

355 

9 

27 

27/28 

+0.10 

2.9 

2.2 

+12.3 

♦38.3 

11.1 

25.4 

30140 

355 

15 

357 

11 

23 

22/23 

+0.10 

3.0 

4.9 

+  6.9 

+35.3 

5.2 

19.6 

30140 

357 

15 

359 

11 

13 

13/12 

+0.25 

3.8 

4.1 

+35.2 

♦63.7 

20.3 

38.5 

30140 

359 

15 

361 

11 

16 

16/14 

+0.20 

3.8 

-15.0 

+79.8 

8.7 

47.8 

30140 

361 

15 

363 

11 

15 

14/14 

♦0.13 

IB 

-16.9 

+57.7 

11.1 

33.7 

30140 

363 

17 

365 

11 

21 

19/21 

+0.11 

-22.6 

+19.5 

14.2 

14.8 

Naan 

+0.10 

D 

m 

1 

14.2 

25.3 

Standard  deviation 

0.09 

0.5 

0.8 

6.8 

10.7 

OOP  -  drag  cat 1m t Ion 
OIP  -  no-drag  aatlMtloo 
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1.  Th«  average  of  Che  orbit  determination  data  realduals 
improved  0.7  m  (17%). 

2.  The  average  of  the  extrapolated  data  residuals  Improved 
11.1  m  (44%) . 

-15  2 

3.  The  mean  bias  force  (0.10  x  10  Rq/s  ^  was  a'>out  w*iat 

we  expected  (compensating  for  the  30  unit  error  in  the  mean  solar 
index) . 

4.  The  worst  drag-estimation  extrapolated  ephemerls  (30  m, 
days  347  to  349)  was  better  than  three  no-drag-estlmatlon  extrapola¬ 
tions  (33.7,  38.5,  and  47.8  m)  and  about  the  same  as  two  other  no- 
drag-estlmatlon  extrapolations  (29.7  and  28.7  m);  that  is,  five  of 
the  standard  runs  did  as  poorly  as  or  worse  than  the  single  worst 
drag-estimation  extrapolated  case. 

Figure  2  summarizes  the  results  of  drag  estimation  versus 
no-drag  estimation.  The  superiority  of  drag  estimation  is  clearly 
evidant.  The  lower  portion  of  Fig.  2  compares  drag  estimation  to 
no-drag  estimation.  The  shaded  portions  outline  the  areas  where 
drag  estimation  was  superior  to  no-drag  estimation. 

It  is  of  Interest  to  speculate  on  the  source  of  the  drag 
model  errors  (fitted  bias  force)  by  looking  for  a  possible  correla¬ 
tion  with  known  variations  in  the  upper  atmosphere. 

Figure  3  is  an  attempt  to  correlate  the  fitted  bias  force 
(correction  for  drag  errors)  with  the  prevailing  solar  and  magnetic 
indices.  First,  we  can  remove  one  known  source  of  error  that  we 
introduced  (30  units  in  the  mean  solar  index),  which  is  equivalent 
—15  2 

to  the  0.10  x  10  Rq/ a  mean  fitted  bias  force.  Only  four  of  the 

12  Intervals  deviate  significantly  from  the  mean:  days  344  to  347 
and  days  357  to  361.  The  earlier  period  coincides  with  a  magnetic 
storm  of  moderate  Intensity,  which  suggests  that  the  modeled  drag 
force  underestimates  the  effects  of  magnetic  storms  (negative  bias 
force).  There  is  no  obvious  correlation  with  either  the  solar  or 
the  magnetic  index  for  the  latter  period  (days  357  to  361). 

The  results  of  the  new  approach  to  drag  estimation  are  excit¬ 
ing  and  potentially  valuable  to  the  Navy  Navigation  Satellite  Sys¬ 
tem  (NNSS).  We  are  now  in  the  midst  of  the  most  disturbed  part  of 
tha  11  yr,  21st  solar  cycle,  when  the  effects  of  Increased  mean  den- 
sltlas  and  unpredlctabla  variations  in  upper  atmospheric  densities 
become  the  rule  rather  than  tha  exception.  Drag  models  do  well  un¬ 
der  normal  circumstances,  but  they  Incorporate  only  known  variations 
and  are  based  on  historical  (20th  solar  cycle)  rather  than  currant 


Difference 

no-drag  estimation  —  drag  estimation  No-drag  estimation  Drag  estimation 


Along-track  bias  force  (Rq/s2  X  1015) 

Solar  index,  s  Magnetic  index,  AD  (result  of  drag  tracking) 
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Fig.  3  Fitted  along-track  bias  force  and  prevailing  solar/magnetic  conditions  for 
satellite  30140. 
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data*  The  next  step  was  to  saa  how  wall  drag  tracking  would  handle 
a  difficult  parlod  Ilka  that  of  April  1978  whan  Transit  satellites 
suffered  serious  extrapolation  errors  as  a  result  of  drag  model 
errors. 
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6.  THE  APRIL  1978  EXPERIMENT 


The  first  100  days  of  1978  were  marked  by  a  sudden  Increase 
In  the  sun*s  activity.  The  mean  solar  index  almost  doubled  between 
November  1977  and  April  1978.  Daily  solar  readings  reached  190* 
and  several  magnetic  storms  of  moderate  intensity  (Ap  >  50)  were 

observed.  During  the  same  period*  the  extrapolated  ephemerides  of 
the  five  Transit  satellites  indicated  that  the  modeled  drag  force 
was  usually  too  low.  The  problem  did  not  surface  until  early  April 
1978  when  a  combination  of  severe  solar  and  magnetic  disturbances 
resulted  in  a  serious  degradation  in  the  accuracy  of  the  predicted 
(extrapolated)  positions  of  all  five  satellites.  A  20  day  span 
(days  90  to  110,  1978)  was  chosen,  and  orbits  were  determined  for 
all  five  satellites  using  the  drag  estimation  option  in  ODP  and  the 
new  approach  described  earlier.  Results  are  summarized  in  Table  5* 
and  the  extrapolated  along-track  navigation  residuals  are  displayed 
in  Figs.  4  through  8.  Orbits  were  determined  for  three  of  the  sat¬ 
ellites  (30120,  30130,  and  30140)  every  other  day  and  for  the  other 
two  daily.  The  worst  extrapolation  errors  occurred  on  days  100  to 
103.  Disturbed  solar  and  magnetic  conditions  on  day  100  are  prob¬ 
ably  responsible  for  the  large  extrapolation  errors  seen  on  days 
100  to  101.  No  such  correlation  exists  for  days  102  to  103,  yet 
the  along-track  extrapolated  positions  were  the  worst  for  all  five 
satellites  (see  Table  5) •  We  will  postpone  the  discussion  of  this 
point  to  the  last  section  of  the  report. 

Figure  9  is  a  plot  of  the  fitted  bias  force  (modeled  drag 
error)  for  each  of  the  five  satellites  and  the  prevailing  solar  and 
magnetic  conditions  for  the  period.  Note  the  following  points  of 
interest: 

1.  The  variations  in  the  fitted  bias  force  (modeled  drag 
force  error)  are  similar  for  all  five  satellites  (see  the  dashed 
line) • 


2.  There  is  a  strong  suggestion  that  the  fitted  bias  force 
follows  the  daily  variation  in  the  magnetic  index,  Ap.  Peak  bias 

forces  were  fitted  out  on  days  100  to  101.  Day  100  recorded  the 
highest  magnetic  index  (A^  *  45).  An  earlier  magnetic  peak  on 

day  93  (A^  -  27)  correlates  well  with  bias  force  peaks  for  all  but 

satellite  30140.  A  late  peak  on  day  103  (A  -  31)  also  correlates 

P 

well  with  (negative)  Increases  in  the  bias  force  on  day  104.  The 
decline  in  the  magnetic  index  on  the  next  three  days  (104  to  106) 
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Fig.  9  Fitted  along-track  bias  forces  and  prevailing  solar/magnetic  conditions. 
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can  also  be  observed  in  the  (negative)  decline  in  the  fitted  bias 
force  on  days  105  to  108. 

3.  The  maximum  (negative)  fitted  bias  force  for  all  three 
satellites  (days  100  to  101)  coincides  with  the  peak  daily  solar 
index  (s  *  180) . 

We  can  draw  two  inferences  from  the  above  observations: 

1.  The  drag  model  underestimates  the  effects  of  Increased 
magnetic  activity,  and 

2.  The  drag  model  underestimates  the  effects  of  daily  vari¬ 
ations  in  the  solar  activity. 

The  poor  performance  of  the  drag  model  during  this  period 
prompted  NAVASTROGRU  to  implement  a  manual  procedure  of  monitoring 
the  extrapolated  along-track  position  errors  daily  and  periodically 
adjusting  the  mean  drag  force  by  changing  the  "true"  (nominal)  co¬ 
efficient  of  drag  parameter  (Cd  -  3.0)  either  upward  (model  under¬ 
predicting  drag)  or  downward  (over-predicting  drag) . 


-  28  - 


THE  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

laurel.  Maryland 


7.  ADDITIONAL  RESULTS  OF  DRAG  ESTIMATION 


Four  additional  experiments  with  drag  estimation  were  carried 
out  over  the  past  year:  (a)  satellite  30200,  days  264  to  304,  1978; 
(b)  satellite  30120,  day  340,  1978,  to  day  15,  1979;  (c)  satellite 
30140,  day  344,  1978,  to  day  15,  1979;  and  (d)  satellite  30190,  days 
94  to  118,  1979.  The  last  three  experiments  were  carried  out  by 
NAVASTROGRU  personnel  as  part  of  their  ODF  and  drag  estimation  op¬ 
tion  checkout. 

The  results  of  the  first  experiment  (satellite  30200)  are 
illustrated  in  Fig.  10.  We  note  the  following: 

1.  The  mean  of  the  rms  of  the  along-track  errors  (48  h  in¬ 
terval,  the  last  24  h  extrapolated)  was  13.4  ±  8.9  m. 

2.  The  fitted  (negative)  along-track  bias  Increased  steadily 
over  the  40  day  span  (21  September  to  31  October),  coinciding  with 
known  deficiencies  in  the  modeled  semiannual  variation  in  density 
(the  maximum  occurring  in  October)  and  the  absence  of  the  seasonal- 
latitudinal  variation  in  helium  content  from  the  presently  imple¬ 
mented  drag  model  (Ref .  2) . 

3.  We  note  a  reasonable  correlation  between  A  and  the  fitted 

P 

bias  force  (-F^)  during  a  moderate  storm  (days  268  to  272) . 

4.  Peak  fitted  bias  forces  on  days  291  to  292  are  not  pre¬ 
ceded  nor  do  they  coincide  with  a  magnetic  storm  (high  A^'s),  al¬ 
though  they  do  occur  two  days  after  maximum  reported  solar  distur¬ 
bances  (8s  180) . 

5.  Manual  drag  estimation  (changing  C^)  Instituted  at 

NAVASTROGRU  following  the  problems  in  April  1978  closely  paralleled 
the  fitted  along-track  bias  force.  varied  from  2.8  (less  than 

nominal)  at  the  start  of  the  span  to  5.4  (1.8  times  nominal)  at  the 
end  (see  item  1  above) .  The  two  local  peaks  in  fitted  bias  force 
(items  3  and  4  above)  have  their  counterparts  in  Cj.  Although  both 


Ref.  2.  A.  Eisner  and  S.  M.  Ylonoulis,  "Long-Period  Terms 
in  the  Neutral  Density  of  the  Upper  Atmosphere,"  JHU/APL  CP  076, 
Nov  1979. 
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methods  follow  the  long-period  variations  well,  the  fitted  bias 
force  is  more  sensitive  to  low-level  errors  and  picks  up  changes 
earlier  than  the  manual  system. 

The  results  of  the  second  experiment  (satellite  30120)  are 
illustrated  in  Fig.  11.  We  note  the  following: 

1.  The  mean  of  the  rms  of  the  along-track  errors  (48  h  in¬ 
terval,  the  last  24  h  extrapolated)  was  21,3  ±  13.5  m,  8  m  worse 
than  the  results  of  the  previous  experiment.  Six  of  the  41  days 
had  an  rms  greater  than  30  m.  Excluding  those  six,  the  rms  of  the 
remaining  35  days  was  16.7  i  6,7  m.  Each  of  the  six  poor  extrapo¬ 
lations  is  associated  with  sudden  changes  to  the  fitted  along-track 
bias  force.  On  three  of  the  six  days,  the  changes  correlate  with 
the  onset  or  end  of  a  magnetic  storm  (days  350  to  351,  1978,  and 
day  4,  1979).  The  balance,  particularly  days  9  and  12,  1979,  cor¬ 
relates  with  large  changes  in  the  fitted  along-track  bias  force 

—16  2 

(>4  x  io  Rq/s  ^  unknown  origin. 

2.  The  fitted  bias  force  shows  a  strong  long-period  (28 
days)  oscillation.  It  is  primarily  a  resonance  effect  (Ref,  3) 
resulting  from  errors  in  the  27th-order  modeled  geopotential  coeffi¬ 
cient.  The  very  pronounced  variations  in  s  and  associated  modeling 
errors  of  this  effect  are  probably  the  other  contributions  to  the 
periodic  variation  in  Ffi. 

3.  There  is  a  possible  correlation  between  A  and  -F„  during 

P  " 

moderate  magnetic  storms  (days  352  to  353) . 

4.  Manual  drag  estimation  (changing  C^)  follows  the  general 

trend  of  the  fitted  bias  force  but  is  less  detailed  (see  the  first 
experiment,  item  5  above). 

The  results  of  the  third  experiment  (satellite  30140)  are 
Illustrated  in  Fig.  12.  We  note  the  following: 

1.  The  mean  of  the  rms  of  the  along-track  errors  (48  h  in¬ 
terval,  the  last  24  h  extrapolated)  was  13.8  ±  9.1  m,  which  is  more 
in  line  with  the  results  of  the  first  experiment.  Only  two  of  the 
37  days  had  rms  values  exceeding  30  m  (excluding  those  two,  the  mean 
was  12.4  ±  7.0  m). 


Ref.  3. S.  M.  Ylonoulls,  "A  Study  of  the  Resonance  Effects 
Due  to  the  Earth's  Potential  Function,"  J.  Geophys.  Res.,  Vol.  70, 
No.  24,  Dec  1965,  pp.  5991-5996,  and  Vol.  71,  No.  4,  Jan  1966, 
pp.  1289-1291. 
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Fitted  along- track  bias  force  (OOP),  operationally  used  Cj  (OIP),  and  prevailing 
solar /magnetic  conditions  for  satellite  30120. 
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Fitted  along-track  bias  force  (ODP),  operationally  used  (OIP),  and  prevailing 
solar/magnetic  conditions  for  satellite  30140. 
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2.  The  strong  periodic  oscillation  in  the  fitted  along- track 
bias  force  noted  with  satellite  30120  (the  second  experiment)  is 
not  evident  in  this  satellite.  Satellite  30140  is  insensitive  to 
errors  in  the  27th-order  modeled  geopotential  coefficient.  The 
weak  periodic  oscillation  that  is  present  is  probably  associated 
with  the  pronounced  variations  in  s  as  discussed  in  the  second  ex¬ 
periment,  item  2. 

3.  There  is  a  strong  suggestion  for  a  correlation  between 
A  and  -Fg  (see  days  348  to  349  and  352  to  353,  1978,  and  days  3 

to  5,  1979). 

4.  Manual  drag  estimation  (changing  C^)  follows  the  general 

trend  of  the  fitted  bias  force  but,  once  again,  it  lacks  the  de¬ 
tailed  structure  of  the  latter  (see  the  discussion  of  the  first  ex¬ 
periment,  item  5). 

The  set  of  experiments  with  drag  estimation  was  completed 
with  the  final  and  most  recent  data  for  satellite  30190.  One  minor 
revision  was  made  in  the  drag  tracking  procedure.  In  all  the  pre¬ 
vious  experiments,  no  explicit  changes  to  the  satellite  period  were 
made;  that  is,  there  was  no  linear  term  in  the  fitting  space.  We 
found  that  we  could  achieve  better  overall  results  by  including  the 
linear  term  in  the  fitting  space  and  "binding"  changes  to  the  semi- 
major  axis  (linear  term)  to  a  few  centimeters.  Results  are  summa¬ 
rized  in  Table  6  and  illustrated  in  Fig.  13.  We  note  the  following: 

1.  The  mean  of  the  rms  of  the  along-track  error  was 
21.3  ±  23.0  m  when  fitting  an  along-track  force,  compared  to 

16.5  ±  16.2  m  when  using  manual  drag  estimation  (changing  C^).  The 

reason  for  the  relatively  poor  performance  of  the  former  (when  com¬ 
pared  to  the  latter)  will  be  dealt  with  in  the  next  section. 

2.  The  fitted  (negative)  along-track  bias  force  decreased 
steadily  over  the  25  day  span  (4  to  28  April) ,  coinciding  with  the 
known  deficiencies  in  the  modeled  semiannual  variation  in  density 
(the  decline  from  the  April  maximum)  and  the  absence  of  seasonal- 
latitudinal  variations  in  atmospheric  helium  content  from  the  pres¬ 
ently  implemented  drag  model  (Ref.  2). 

3.  A  strong  correlation  is  evident  between  A^  and  -Fg  on 

three  separate  occasions  (days  95,  112,  and  116).  In  all  three  the 
modeled  drag  force  is  too  small  and  is  supplemented  by  sudden  in¬ 
creases  in  the  fitted  along-track  bias  force  (days  95,  113,  and  116). 


Table  6 
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Satellite  30190  in  1979. 


Day 

Solar 

index 

Magnetic 

Index 

A 

P 

No-drag  tracking 

Drag  tracking 

3 

month 

mean 

Dally 

ca 

Along- track  rma  (a) 

-Fb  *  1015 
(C4  -  3.0) 

Along-track 

»  <m) 

Prediction 

Track 

Prediction 

Track 

94 

198 

200 

4.5 

18.2 

4.7 

■ 

3.2 

95 

198 

183 

4.5 

84.0 

17.5 

wBM 

1 

5.5 

96 

198 

179 

13 

4.7 

10.1 

6.0 

§■ 

3.8 

97 

197 

171 

9 

4.7 

40.0 

13.5 

0.8 

21 

1 

3.2 

96 

197 

mm 

8 

4.0 

16.7 

7.3 

:  V-B 

19 

9 

3.2 

99 

197 

■a 

8 

3.8 

10.2 

3.5 

1 

5 

2 

3.5 

100 

197 

170 

10 

3.8 

4.8 

3.7 

1 

11 

2 

3.1 

101 

197 

179 

5 

4.2 

9.4 

3.9 

I 

17 

3 

2.9 

102 

197 

170 

9 

4.2 

1  ■ 

3.9 

|  '1 

Eli 

3.3 

103 

197 

172 

9 

4.2 

2.8 

9 . 

2.1 

104 

197 

170 

11 

4.2 

4.5 

0.8 

2 

3.3 

105 

196 

167 

15 

3.8 

1 

5.2 

0.4 

15 

9 

2.8 

106 

196 

172 

12 

3.8 

2.8 

0.7 

11 

5 

2.9 

107 

196 

170 

12 

3.7 

9.4 

4.4 

0.5 

20. 

8 

2.5 

108 

196 

157 

6 

3.7 

6.1 

3.1 

0.5 

5. 

2.5 

109 

195 

152 

4 

3.5 

6.8 

4.1 

0.2 

■ 

1 

2.7 

110 

195 

151 

1 

3.1 

26.5 

4.4 

0.2 

3.0 

111 

194 

165 

12 

3.1 

19.7 

5.3 

0.3 

2.8 

112 

193 

159 

3.3 

5.5 

3.2 

0.3 

5. 

6 

2.9 

113 

192 

162 

18 

3.6 

- 

- 

0.8 

34 

1 

2.8 

114 

192 

163 

14 

3.6 

8.3 

7.1 

0.2 

ULU 

3.6 

115 

191 

174 

r*9i 

3.3 

11.3 

4.3 

0.2 

8 

7 

3.6 

116 

191 

178 

11 

3.7 

18.2 

10.0 

1.2 

42 

0 

3.7 

117 

191 

185 

22 

3.6 

13.7 

11.4 

-0.3 

Lzl 

i\ 

3.4 

118 

191 

182 

31 

3.1 

16.1 

5.9 

-j 

-0.2 

4 

4 

3.4 

Naan* 

16.5  ±  16.2 

5.7  ±  3.8 

21.3 

23.0 

3.2  t  0.6 

♦Exclude*  day  113 
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Day  (1979) 


Fig.  13  Fitted  along-track  bias  force  (ODP),  operationally  used  C<j  (OIP),  and  prevailing 
solar/magnetic  conditions  for  satellite  30190. 
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4.  Manual  drag  estimation  (changing  C^)  follows  the  general 

trend  of  the  fitted  bias  force  but  lacks  the  detailed  structure  of 
the  latter  (see  the  discussion  of  the  first  experiment.  Item  5). 
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8.  REFINING  THE  DRAG  ESTIMATION  TECHNIQUE 


The  results  of  the  last  experiment  described  In  the  previous 
section  were  surprising.  We  expected  that  the  computer  drag  esti¬ 
mation  (fitting  a  bias  force)  would  do  better  than  its  manual  coun¬ 
terpart  of  changing  C^.  Instead  we  found  that  the  opposite  was 

true.  Let  us  examine  the  assumptions  that  form  the  basis  for  com¬ 
puter  drag  estimation. 

The  first  assumption  Is  that  removing  the  quadratic  growth 
in  the  along-track  position  error  will  result  In  a  more  precise  set 
of  initial  orbit  elements  that,  when  integrated,  will  give  a  better 
fit  to  the  data  (along-track  errors). 

The  second  assumption  is  that  errors  In  the  drag  model  are 
slowly  varying  and  that  the  best  guess  at  tomorrow's  force  correc¬ 
tion  Is  to  use  today's  fitted  bias  force. 

The  first  assumption  Is  clearly  correct.  The  mean  of  the  rms 
of  the  along-track  errors  In  the  tracking  span  of  the  last  experi¬ 
ment  (Table  6)  was  3.2  ±  0.6  m  as  compared  to  5.7  ±  3.8  m,  a  43% 
reduction  In  favor  of  computer  drag  estimation.  The  second  assump¬ 
tion  is  only  partially  correct;  that  is,  although  it  is  generally 
true  that  drag  model  errors  are  overall  slowly  varying  (see  Figs.  9 
through  13),  it  is  not  always  true  that  the  best  guess  at  tomorrow's 
errors  is  today's  fitted  error.  The  two  worst  predictions  in  the 
last  experiment  (satellite  30190)  on  days  96  and  117  (shaded  areas 
in  Fig.  13)  are  a  result  of  using  the  large  (transient)  fitted 
along-track  bias  force  on  days  95  and  116  to  predict  (extrapolate) 
the  satellite's  along-track  position  on  days  96  and  117,  respec¬ 
tively.  There  are  several  other  instances,  although  less  dramatic, 
in  the  earlier  experiments  described  above.  One  example  was  satel¬ 
lite  30120,  day  350,  1978,  and  the  extrapolated  position  errors  on 
the  following  day  (shaded  areas  in  Fig.  11).  Another  example  was 
the  large  though  transient  along-track  bias  forces  (Fig.  9)  fitted 
for  days  100  to  101,  1978,  associated  with  solar  and  magnetic  dis¬ 
turbances  that,  when  used  in  the  prediction  span  (days  102  to  103), 
resulted  in  large  along-track  position  errors  (shaded  areas  in  Figs. 
4  to  8). 

In  every  case  we  are  faced  with  a  large  change  (greater  than 
"16  2 

4  x  10  R  /s  )  in  the  fitted  along-track  bias  force  from  one  track 
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to  the  next.  The  fitted  force  is  correct  for  the  track  (orbit  de¬ 
termination  Interval) ,  but  it  is  likely  to  result  in  large  along- 
track  extrapolated  position  errors  since  the  source  of  the  change 
was  almost  invariably  transient  in  nature  (a  magnetic  storm)  and, 
as  a  rule,  does  not  persist  beyond  the  tracking  span. 

A  safer  strategy  for  dealing  with  rapid  changes  in  the  fitted 
along-track  bias  force  is  to  replace  force  with  a  weighted  mean  of 
the  preceding  two  to  five  fitted  forces  (weighting  the  most  recent 
most  heavily,  etc.).  Using  such  a  strategy,  days  96  and  117  (sat¬ 
ellite  30190,  1979,  Table  6)  were  repeated  with  the  fitted  forces 

of  27  x  10  R^/s2  replaced  with  19  x  10  ^  R^/s2  *-n  t*le  former  and 

12  x  io”*^  Rq/s2  replaced  with  7  x  10  ^  RQ/s2  in  t*ie  latter.  The 

rms  of  the  prediction  Improved  from  86.1  to  23  m  in  the  former  and 
from  74.9  to  22  m  in  the  latter.  The  new  mean  of  the  rms  of  the 
along-track  errors  was  16.4  ±  13.5  m,  a  significant  Improvement  of 
5  m  in  the  mean  and  10  m  in  the  standard  deviation.  In  practice, 
the  strategy  described  is  simple  to  implement.  It  is  only  necessary 
to  keep  a  dally  log  of  the  fitted  along-track  bias  force  for  each 
satellite.  As  a  new  orbit  is  determined,  the  latest  fitted  force 
is  checked  against  the  log.  Whenever  the  fitted  force  (new  force) 
differs  from  the  input  force  (previous  run  fitted  force)  by  more 
—16  2 

than  4  x  io  Rq/s  (in  magnitude),  it  is  replaced  with  a  weighted 

mean  of  the  previous  three  to  four  determined  forces  (from  the  log), 
and  a  new  extrapolated  ephemeris  is  generated  for  subsequent  injec¬ 
tion  into  the  satellite  memory. 
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Appendix  A 

DRAG  ESTIMATION  THEORY* 


We  begin  with  a  set  of  so-called  HLC  equations  of  satellite 
motion  that  were  derived  by  W.  H.  Guier: 


H"  -  (3+10  e  cos  M)H  -  (2+4  e  cos  M)L*  +  (2  c  sin  M)L  * 

n 


L"  +  (2+4  e  cos  M)Hf  -  (2  e  sin  M)H  -  (e  cos  M)L  *=  ~ 

n 


C"  +  (1+3  e  cos  M)C  =  , 

n 


(A-l) 


where 


H, L,C  are  the  range,  along-track,  and  cross-track  position 

errors,  and  H",L" ,C"  are  the  respective  associated 
perturbation  accelerations , 

F  ,F  ,F  are  the  perturbing  accelerations  in  the  H,L,C 
h  L  L  directions, 

e  is  the  satellite  eccentricity, 

n  is  the  satellite  mean  motion,  and 

M  is  the  satellite  mean  anomaly* 

We  make  the  following  assumptions: 

I,  The  orbits  are  circular,  i.e.,  e  *  0;  and 

2*  The  force  of  drag  is  acting  along-track  only,  i.e*, 


^Developed  by  S.  M.  Yionoulis 
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Using  these  assumptions,  Eq,  A-l  becomes 
H"  -  3H  -  2L1  «  0 

fl 

L"  +  2H'  =  -j 
n 

C"  +  C  =  0  .  (A-2) 


We  now  Integrate  the  second  equation  in  Eq.  A-2, 


L'  +  2H  ■  — j  , 


(A-3) 


and,  substituting  into  the  first  equation  in  Eq.  A-2,  we  obtain 


H"  -  3H  -  2j  —j  M  -  2Hl 

V" 

fl 

H'  +  -2  M  =  0  . 

n 


(A-4) 


The  solution  of  Eq.  A-4  is 


fl 

H  -  2  -j  M  . 
n 


(A-5) 


Substitute  Eq.  A-5  into  Eq.  A-3, 
F  F 

L'  +  4-~M-~M  , 

n  n 

and  by  integrating  we  obtain 

L  -  -3/2  -j  M2 
n 


H  -  2  -j  M  . 
n 


(A-6) 
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The  next  step  Is  to  arrive  at  an  expression  relating  to 

one  of  the  drag  parameters  (C^) .  Is  the  error  In  the  modeled 

drag  force.  We  do  not  propose  abandoning  the  drag  model  but  rather 
supplementing  It  with  drag  estimation. 

We  can  write  the  following  relation  for  the  error  In  the 
modeled  drag  acceleration,  : 


Fl  -  1/2  A/M  (p>  V^Cj-Cj  ) 


(A-7) 


where 


A/M  is  the  area-to-mass  ratio  for  the  satellite, 

(p)  Is  the  average  density  over  the  tracking  span, 

V  is  the  average  satellite  velocity  over  the  tracking 
span, 

C.  Is  the  input  coefficient  of  drag,  and 
d0 

Cj  is  the  adjusted  (tracked)  coefficient  of  drag. 


The  actual  parameter  used  in  the  fitting  process  (a^)  is 
related  to  F^  as  follows: 


a6 


(A-8) 


and  the  correction  to  using  Eqs.  A-7  and  A-8  is 


-4/3 


2 

n  a6  § 
A/M  <p>  V2 


The  concept  of  drag  estimation  was  later  generalized  to 
fitting  out  an  along-track  bias  force  Ffi  and  Integrating  the  force 

explicitly  rather  than  implicitly  via  the  change  In  C^.  Note  that 
changing  Is  not  the  same  as  using  Ffi.  The  latter  (Ffi)  Is  a  con¬ 
stant  force,  whereas  the  former  Is  a  varying  (following  the  vari¬ 
ations  in  p)  force  whose  mean  approximates  Fg. 
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